>e type a plus sign (+} in 



n' 1 ITII ITV 

UTILITY 


Attorney Docket No. 1 1 5-CT-527 1 


IE<" PATPMT ADPI ir*ATir>M 
g= r M I CIn I ArrLluA 1 IUiM 


F;rsr Named ln\i 


entor or Application Identifier | Min Xie et al. 


TRANSMITTAL 


Title 


METHOD AND APPARATUS FOR FAST NATURAL LOGfX 
CALCULATION 


jjpg'jjj - new nonprovisional applications under 37 CFR 1.53(b)) 


Express Mail Label No |EL31 9728805US 



ADDRESS TO: 



Washington, DC 20231 



^1 



jX I Fee Transmittal Form (e.g., PTO/SB/17) 

(Submit an original, and a duplicate for fee processing) 

Specification [Total Page 

(Preferred arrangement set forth below) 

- Descriptive title of the Invention 

- Cross References to Related Applications 

- Statement Regarding Fed sponsored R&D 

- Reference to Microfiche Appendix 

- Background of the Invention 

- Brief Summary of the Invention 

- Brief Description of the Drawings (if filed) 

- Detaiied Description 

- Claim(s) 

- Abstract of the Disclosure 
3. H Drawing(s) (35 USC 11 3) [Total Sheets 

Oath or Declaration [Total Pages 

a. | | Newly executed (original or copy) 

b. | | Copy from a prior application (37 CFR 1 63(d)) 

(for continuation/divisional with Box 17 completed) 

. . [Note Box 5 below] 

I. | I DELETION OF INVENTOR(S) 

Signed statement attached deleting inventor(s) named 
in the prior application, see 37 CFR 1 .63(d)(2) and 
1.33(b). 

j | Incorporation by Reference (useable if Box 4b is checked) 

The entire disclosure of the prior application, from which a copy of the 
oath or declaration is supplied under Box 4b, is considered as being part 
of the disclosure of the accompanying application and is hereby 
incorporated by reference therein. 



□ 



Microfiche Computer Program (Appendix) 



Computer Readable Copy 

Paper Copy (identical to computer copy) 

Statement verifying identity of above copies 



[O 

n 



ACCOMPANYING APPLICATION PARTS 



HAssignmei 
37 CFR 3. 



Assignment Papers (cover sheet & documents)) 

| | Power of Attorney 



| [ English Translation Document (if applicable) 

Information Disclosure Copies of IDS 

1 1 Statement (I DS)/PTO-1 449 ' 1 Citations 

Preliminary Amendment 



□ s 



Small Entity Statement filed in prior application, 

Statement(s) Status still proper and desired 

(PTO/SB/09-12) 

5. Certified Copy of Priority Document(s) 
(If foreign priority is claimed) 

Other: Express Mail Certificate; Declaration 
purposes of identification of inventors only 



Power of Attorney for 



•WO TE FOR ITEMS 1 &U 



APPLICATION IS RELIED UPON (37 C F i 



17. Ifa CONTINUING APPLICATION, check appropriate box and supply the requisite information: 

_| Continuation \~\ Divisional I I Continuation-m-partfCIP) 

ifiorapplicaton information L EX&mmer brdip/Art Unit 



of pnor application No. 



18. CORRESPONDENCE ADDRESS 



ICustomer Number or Bar Code Label 



(Insert Customer No ar Attach ba 



□ 



Correspondence address below 



Armstrong Teasdale LLP 



One Metropolitan Square 



CITY 


St. Louis 


| STATE 


MO 


ZIP CODE 


163102-2740 




COUNTRY 


U.S.A. 


I TELEPHONE 


(314) 621-5070 


FAX 1(314)621-5065 




Name (Print/type) 


Alan L. Cassel 


Registration No. (Attorney/Agent) 


35,842 






Signature 




| Date 







!!f 0 «™ ,atement: ™ S f T is est i m 4 at i d to f take 02 hours to complete. Time will vary depending upon the needs of the individual case. Any comments on the 
amount of time you are required to complete this form should be sent to the Chief Information Office, Patent and trademark Office. Washinaton DC 20231 DO NOT 



PATENT 
15-CT-5271 



CERTIFICATE OF MAILING BY EXPRESS MAIL TO 
THE COMMISSIONER OF PATENTS AND TRADEMARKS 



Express Mail mailing label number: EL319728805 US 
Date of Mailing: February 18, 2000 



I certify that the attached continuation patent application of MIN XIE AND 
SHARON X. WANG for METHOD AND APPARATUS FOR FAST NATURAL 
LOG(X) CALCULATION (Attorney Docket No.l5-CT-5271) including: 

• Certificate of Mailing Via Express Mail (1 pg.) 

• Utility Patent Application Transmittal (1 pg.) 

• Fee Transmittal (in duplicate) (1 pg.) 

• Declaration and Power of Attorney (2 pgs.)(for purposes of identification 
only) 

• Eight (8) pages of specification; eight (8) pages of claims; one (1) page of 
abstract 

• Two (2) sheets of drawings 

• Return post card 

is being deposited with the United States Postal Service "Express Mail Post Office to 
Addressee" service under 37 C.F.R. §1.10 on the date indicated above in an envelope 
addressed to the Assistant Commissioner for Patents, Box PATENT APPLICATION, 
Washington, D.C. 20231. 




Reg. No.3 5,842 

Armstrong Teasdale LLP 

One Metropolitan Square, Suite 2600 

St. Louis, MO 63102 

(314) 621-5070 



Express Mail No. EL319728805US 



METHOD AND APPARATUS FOR FAST 
NATURAL LOG(X) CALCULATION 

BACKGROUND OF THE INVENTION 

The invention relates generally to methods and apparatus for 
computing a computationally intensive algorithm, and more specifically to a method 
and apparatus for computing log(x), or equivalently, -log(x), in a manner that is 
particularly useful for computed tomographic image processing and other 
5 applications. 

In at least one known computed tomography (CT) imaging system 
configuration, an x-ray source projects a fan-shaped beam which is collimated to lie 
within an X-Y plane of a Cartesian coordinate system and generally referred to as the 
"imaging plane." The x-ray beam passes through the object being imaged, such as a 

10 patient. The beam, after being attenuated by the object, impinges upon an array of 
radiation detectors. The intensity of the attenuated beam radiation received at the 
detector array is dependent upon the attenuation of the x-ray beam by the object. Each 
detector element of the array produces a separate electrical signal that is a 
measurement of the beam attenuation at the detector location. The attenuation 

15 measurements from all the detectors are acquired separately to produce a transmission 
profile. 

In known third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged so that the angle at which the x-ray beam intersects the object constantly 

20 changes. A group of x-ray attenuation measurements, i.e., projection data, from the 
detector array at one gantry angle is referred to as a "view." A "scan" of the object 
comprises a set of views made at a different gantry angles, or view angles, during one 
revolution of the x-ray source and detector, hi an axial scan, the projection data is 
processed to construct an image that corresponds to a two dimensional slice taken 

25 through the object. One method for reconstructing an image from a set of projection 



-1- 



15-CT-5271 



data is referred to in the art as the filtered back projection technique. This process 
converts the attenuation measurements from a scan into intergers called "CT numbers" 
or "Hounsfield units," which are used to control the brightness of a corresponding 
pixel on a cathode ray tube display. 

5 The negative natural logarithm function -log(x) is an important yet 

computationally intensive algorithm in computed tomographic (CT) image processing. 
In known systems, a 5th order polynomial is used to approximate the function. 
However, this polynomial still consumes more than 20% of the total image processing 
time and generates a relatively large approximation error and error standard deviation. 

10 A positive floating point number x can be represented by an expression 

written as: 

x = mx2 e (1) 

where m (I < m <2) is a. mantissa and e is a binary exponent. 

Using equation (1), -log(x) can be written as: 

15 y = - log(x) = - log(m) -ex log(2) (2) 

The following equation uses a finite order polynomial to approximate 
log(m) in a region 1 < m < 2 . Generally speaking, the higher the order of the 
polynomial, the better the approximation will be, but the computational load is in 
proportion to the order of the polynomial. For example, a 5th order polynomial 
20 presently used is written as: 

log(m) « (a 0 +a l m +a 2 m 2 +a 3 m 3 +a 4 m 4 +a 5 m 5 ) (3a) 

or as: 

y = -log(x) « -(a 0 + m{a x + m(a 2 + m{a^ + m(a 4 + a 5 m)))) + e x a 6 ) (3b) 
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In equation (3b), a 0 ~ a s are precalculated constants. To compute 
-log(x), six additions and six multiplications are required, plus mantissa and exponent 
extractions. 

To process images more efficiently and accurately, it would be 
5 desirable to provide methods and apparatus to reduce the complexity of the 
approximation used to calculate -log(x) while achieving numerical accuracy consistent 
with IEEE (Institute of Electrical and Electronic Engineers) floating-point precision. 

BRIEF SUMMARY OF THE INVENTION 

There is therefore provided, in one embodiment of the present 
invention, a method for computing a natural logarithm function that includes steps of: 

10 partitioning a mantissa region between 1 and 2 into N equally spaced sub-regions; 
precomputing centerpoints a t of each of the N equally spaced sub-regions, where 
i = 0,...,N -1; selecting N sufficiently large so that, for each sub-region, a first degree 
polynomial in m computes log(m) to within a preselected degree of accuracy for any m 
within the sub-region, where m is a mantissa of a binary floating point representation 

15 of a number; and computing a value of \og(x) for a binary floating point representation 
of a particular number x stored in a memory of a computing device utilizing the first 
degree polynomial in m. 

It will be seen that this embodiment and others described herein reduce 
the complexity of approximations used to calculate natural logarithms while achieving 
20 numerical accuracy consistent with IEEE floating point precision. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a pictorial view of a CT imaging system. 

Figure 2 is a block schematic diagram of the system illustrated in 

Figure 1. 
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Figure 3 is a representation of a number stored in IEEE single- 
precision binary floating point format, partitioned as in one embodiment of the 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Referring to Figures 1 and 2, a computed tomography (CT) imaging 
5 system 10 is shown as including a gantry 12 representative of a "third generation" CT 
scanner. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a 
detector array 18 on the opposite side of gantry 12. Detector array 18 is formed by 
detector elements 20 which together sense the projected x-rays that pass through an 
object 22, for example a medical patient. Detector array 18 may be fabricated in a 
10 single slice or multi-slice configuration. Each detector element 20 produces an 
electrical signal that represents the intensity of an impinging x-ray beam and hence the 
attenuation of the beam as it passes through patient 22. During a scan to acquire x-ray 
projection data, gantry 12 and the components mounted thereon rotate about a center 
of rotation 24. 

15 Rotation of gantry 12 and the operation of x-ray source 14 are 

governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 

20 analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

25 Computer 36 also receives commands and scanning parameters from 

an operator via console 40 that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the reconstructed image and other data from 
computer 36. The operator supplied commands and parameters are used by computer 
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36 to provide control signals and information to DAS 32, x-ray controller 28 and 
gantry motor controller 30. In addition, computer 36 operates a table motor controller 
44 which controls a motorized table 46 to position patient 22 in gantry 12. 
Particularly, table 46 moves portions of patient 22 through gantry opening 48. 

A negative natural logarithm function -log(/») is used by image 
reconstructor 34 to generate images. In one embodiment of the invention, the function 
log(m) of equation (3a) above is written as: 



iogO)'L a , , iog(»"i =n / ^ 2 

log(m) = log(a) + * V i! Mm - c (m - a) + {m-af +■■■ 

-\ ! (m — a) 



(4a) 



log(m) « log(a) + - + • • ■ (4b) 

a 2a 



where a is a known reference point. The error of the above function is 



log(m) ( " +) , 

error < ■ (m - a) 

Nl 



(5) 



Because (m-a)<l, there are two ways to minimize the error. One way 
is to increase the order of the approximation, and the other is to minimize the distance 
from m to a. Because mantissa m is between 1 and 2, in one embodiment of the 
present invention, the region between 1 and 2 is partitioned into N equally spaced sub- 
regions. Centers of each of the sub-regions are precomputed and used as reference 
points in equations (4a) and (4b). By partitioning into a sufficiently large number of 
sub-regions, a low order polynomial function produces sufficient accuracy for CT 
imaging purposes. In particular, by selecting a sufficiently large number of sub- 
regions, for any m within any particular sub-region, log(m) is computed by a first- 
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degree polynomial to within a preselected degree of accuracy within that sub-region. 
For example, computer 36 uses the first degree polynomial in m to compute values of 
log(x) for binary floating point representations of particular numbers x stored in its 
memory. 

5 A log(/ra) approximation that is based on a first order polynomial with a 

set of precalculated reference points is written as follows: 

log(m)«log(a.) + (OT ~ a ' ) ; i = 0,...,N-V, l<a t <2 (6) 
a i 

where a t is a closest reference point to a given mantissa m. 

Rather than compute a sub-region index using i = round ((m -l)xN), 

10 which would require six operations, one embodiment of the present invention reduces 
computation load as follows. A partitioning algorithm divides the mantissa of a 
binary floating point number in memory into two sub-regions. The sub-regions have 
index i and Ax, where Ax is a distance from mantissa m to reference point a { . Indices i 
and Ax are directly extracted from an IEEE floating-point number stored in a computer 

15 system, thereby reducing computation time and improving accuracy. In one 
embodiment, mantissa partitioning occurs as illustrated in Figure 3, in which index i 
ranges from 0 to 127 and each region represents information extracted from the datum 
shown in Figure 3. More particularly, in a single precision IEEE floating point 
number, b 31 represents a sign bit, b 30 the most significant bit of exponent e, b 2 3 the 

20 least significant bit of exponent e, b 22 the most significant bit of mantissa m, and b 0 
the least significant bit of mantissa m. (If it is desired to use a different designation 
for the numbering of bits b, those skilled in the art can make the appropriate changes 
required in the description for notational consistency.) In this single precision 
embodiment, exponent e is extracted directly from bits b 30 to b 23 ; region i is extracted 

25 directly from bits b 22 to b J6 ; and Ax (a distance from mantissa m to reference point ai) 
is extracted directly from bits b 15 to b 0 . 
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Using the extraction illustrated in Figure 3, a maximum error of 
equation (6) in each sub-region is estimated by an expression written as: 

error <JLj-±S) ■ z=0,...,iV-l; \<a t <2 (7a) 
2a, 2 [2NJ 

From equation (7a), it is seen that the error of the first order 
5 approximation is always positive, so that the error is biased. To minimize the 
maximum error, in one embodiment the mean error of equation (7a) is subtracted from 
equation (6). Thus, the unbiased error is written as: 

\error\<{— [ — ] ; i = 0,...,N - 1; l<a,<2 (7b) 
{4Na, ) 

Subtracting equation (7b) from equation (6) results in an unbiased first 
10 order polynomial function for -log(x) written as: 

y = - log(x) ~ b f + c ; - Ax + e x log(2) (8) 
for * = 0,...,iV-l 

c t =-l/a i 

where a = 1 + 1 + Q-5 and Ax is a distance from mantissa m to 

N 

15 reference point a,-. The value Ax is extracted directly from an IEEE floating point 
datum. In one embodiment, log(2) and the b { and a are pre-calculated and saved in a 
look-up table at initialization time. In one embodiment, the values b, are determined 
from precomputed values of log(^ 2 ). For purposes of comparison, computation of 
equation (8) requires only 1/3 as much time as is required to calculate the 5th order 

20 approximation of equation (3b). 
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Taking into account the relationship of m and Ax, an approximation to 
log(x) (or equivalently -log(x)) is computed using a polynomial of first degree in m 
and a precomputed value of log(<3,)- 

In one embodiment of the present invention, image reconstructor 34 is 
5 configured with software or firmware to compute logarithms using one or more of the 
methods of the present invention, when CT imaging system 10 images an object from 
acquired projection data. 

From the preceding description of various embodiments of the present 
invention, it is evident that the complexity of the approximation used to calculate 

10 -log(;t) is reduced, while numerical accuracy consistent with IEEE (Institute of 
Electrical and Electronic Engineers) floating-point precision is maintained. Thus, 
images processed by CT imaging system 10 are processed more efficiently, and 
without loss of detail. Although particular embodiments of the invention have been 
described and illustrated in detail, it is to be clearly understood that the same is 

15 intended by way of illustration and example only and is not to be taken by way of 
limitation. For example, embodiments of the improved computation for log(x) and 
-log(x) can be incorporated into any computational system requiring increased 
efficiency while maintaining computational accuracy. In addition, the present 
invention is suitable for use with floating point numbers having greater or lesser 

20 precision than those discussed in detail in this description. The modifications 
necessary to accommodate such different precisions will be apparent to those skilled 
in the art, once the invention described herein is thoroughly understood. Accordingly, 
the spirit and scope of the invention are to be limited only by the terms of the 
appended claims and legal equivalents. 
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WHAT IS CLAIMED IS: 

1. A method for computing a natural logarithm function comprising 

the steps of: 

partitioning a mantissa region between 1 and 2 into N equally spaced 

sub -regions; 

precomputing centerpoints a t of each of the N equally spaced sub- 
regions, where i - 0, . . . , N - 1 ; 

selecting N sufficiently large so that, for each sub-region, a first degree 
polynomial in m computes log(/w) to within a preselected degree of accuracy for any m 
within the sub-region, where m is a binary mantissa of a binary floating point 
representation of a number; and 

computing a value of \og{x) for a binary floating point representation of 
a particular number x stored in a memory of a computing device utilizing the first 
degree polynomial in m. 

2. A method in accordance with Claim 1 wherein the particular 
number x has a binary exponent e in addition to the binary mantissa m; 

and further wherein computing a value of log(x) for the binary floating 
point representation of the particular number x comprises the steps of: 

partitioning a mantissa m of a binary representation of x in a memory, 
the representation of x including a binary exponent e and the binary mantissa m, 
wherein a first, most significant part of the partition corresponds to a region i and a 
second, less significant part of the partition corresponds to a region Ax, where Ax is a 

t i + 0.5 

distance from mantissa m to reference point a i = 1 + — — — ; and 

computing an approximation to log(x), using a polynomial of first 
degree in m and a precomputed value of log(a r ). 
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3. A method in accordance with Claim 2 wherein computing the 
approximation to log(jc) comprises the step of computing an approximation written as: 

(m -a t ) 
log<» * log(a ; ) + ; 

where a t is a closest reference point to the binary mantissa m of the 

5 number x; and 

\<a i <2. 

4. A method in accordance with Claim 2 wherein computing an 
approximation to log(x) comprises the step of computing an approximation written as: 

y = -log(x) « b i + c l Ax + e x log(2) 

10 for i = 0,...,N-\ 

where: 

b =-log(a I -) + f— — I — ;and 

&K l) {4a { N ) { 2N) ai 

c t =-\la r 

5. A method in accordance with Claim 4 further comprising the steps 
of precomputing a value for log(2), and, for each *, precomputing each value of b, and 

15 a. 

6. A method in accordance with Claim 5 further comprising the step of 
storing the precomputed values of b, and c, in a look-up table. 

7. A method in accordance with Claim 2 wherein the number x is 
represented by a 32-bit representation having a sign bit, an 8-bit exponent, and a 

20 23-bit binary mantissa m having bits b 22 to b 0 in order of significance with 622 heing a 
bit of greatest significance; and the step of partitioning the mantissa m comprises the 

-10- 



15-CT-5271 



step of selecting a first group of bits b 22 through b\$ as index / and bits b\ 5 through b 0 
as Ax. 

8. A method in accordance with Claim 1 utilized in a computed 
tomography (CT) scanner for generating an image of an object from acquired 

5 proj ection data of the obj ect. 

9. A method in accordance with Claim 8 wherein said natural 
logarithm is used in an image reconstructor to generate the image of the object. 

10. A method in accordance with Claim 8 wherein the particular 
number x has a binary exponent e in addition to the binary mantissa m; 

1 0 and further wherein computing a value of log(x) for the binary floating 

point representation of the particular number x comprises the steps of: 

partitioning a mantissa m of a binary representation of x in a memory, 
the representation of x including a binary exponent e and the binary mantissa m, 
wherein a first, most significant part of the partition corresponds to a region i and a 
15 second, less significant part of the partition corresponds to a region Ax, where Ax is a 

1 r + 0.5 

distance from mantissa m to reference point a j = 1 H — — ; and 

computing an approximation to log(x), using a polynomial of first 
degree in m and a precomputed value of log(a ; ). 

11. A method in accordance with Claim 10 wherein computing the 
20 approximation to log(x) comprises the step of computing an approximation written as: 

log(w) « log(a, ) + — — — ; 

a i 

where a t is a closest reference point to the mantissa m; and 
\<a l <2. 
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12. A method in accordance with Claim 10 wherein computing an 
approximation to log(x) comprises the step of computing an approximation written as: 

y = - log(x) ~ b { + c i Ax + e x log(2) 

for i = 0,...,N-l 



-log(a) + [ — - — I -f 1 + — | — ;and 
^ l} {4a ,N ) { 2N)a ; 

-\!a ; . 



13. A method in accordance with Claim 12 further comprising the 
steps of precomputing a value for log(2), and, for each i, precomputing each value of 
bi and ci. 

14. A method in accordance with Claim 13 further comprising the step 
of storing the precomputed values of bi and c t in a look-up table. 

15. A computing device comprising a memory in which binary 
floating point representations of particular numbers are stored, said device being 
configured to: 

partition a mantissa region between 1 and 2 into N equally spaced sub- 
regions; 

precompute centerpoints a, of each of the N equally spaced sub- 
regions, where i=0,...,N-l, wherein N is sufficiently large so that, within each sub- 
region, a first degree polynomial in m computes log(m) to within a preselected degree 
of accuracy for any m within the sub-region, where m is a binary mantissa of a binary 
floating point representation of a number; and 

compute a value of log(x) for a binary floating point representation of a 
particular number x stored in said memory utilizing the first degree polynomial in m. 
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16. A computing device in accordance with Claim 15 wherein the 
particular number x has a binary exponent e in addition to the binary mantissa m; 

and wherein said device being configured to compute a value of log(x) 
for the binary floating point representation of the particular number x comprises said 
device being configured to: 

partition a mantissa to of a binary representation of x in a memory of 
said device, the representation of x including a binary exponent e and the binary 
mantissa m, wherein a first, most significant part of the partition corresponds to a 
region i and a second, less significant part of the partition corresponds to a region Ax, 

1 / + 0.5 , 

where Ax is a distance from mantissa m to reference point a,. =1+ ; and 

compute an approximation to log(x), using a polynomial of first degree 
in m and a precomputed value of log(a i ). 

17. A computing device in accordance with Claim 16 wherein said 
device being configured to compute the approximation to log(x) comprises said device 
being configured to compute an approximation written as: 

(m — a 1 ) 
log(m) « logO,. ) + ; 

where a,- is a closest reference point to the binary mantissa m of the 

number x; and 

1 < a. < 2 . 

18. A computing device in accordance with Claim 16 wherein said 
device being configured to compute an approximation to log(x) comprises said device 
being configured to compute an approximation written as: 

y = -log(jc) « b ; + c, Ax + e x log(2) 
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for i = 0,...,JV-l 
where: 

6 -=- log(a ' )+ ^) "( 1+ ^)i ;md 

c,. =-l/fl f . 

19. A computing device in accordance with Claim 18 further 
configured to precompute a value for log(2), and, for each i, to precompute each value 
of bt and q. 

20. A computing device in accordance with Claim 19 further 
configured to store the precomputed values of b { and c, in a look-up table. 

21. A computing device in accordance with Claim 16 wherein the 
number x is represented by a 32-bit representation having a sign bit, an 8-bit exponent, 
and a 23 -bit binary mantissa m having bits b 22 to b 0 in order of significance with b 2 2 
being a bit of greatest significance; and wherein said device being configured to 
partition the mantissa m comprises said device being configured to select a first group 
of bits b 2 2 through b\e as index i and bits b\ 5 through b 0 as Ax. 

22. A computing device in accordance with Claim 15 in a computed 
tomography (CT) scanner and utilized by said CT scanner for calculating logarithms 
when said CT scanner generates an image of an object from acquired projection data 
of the object. 

23. A computing device in accordance with Claim 22 wherein said CT 
scanner utilizes said computing device to calculate natural logarithm in an image 
reconstructor to generate the image of the object. 

24. A computing device in accordance with Claim 22 wherein the 
particular number x is stored with a binary exponent e in addition to the binary 
mantissa m; 
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and further wherein said device being configured to compute a value of 
log(jc) for the binary floating point representation of the particular number x comprises 
said device being configured to: 

partition a mantissa m of a binary representation of x in a memory, the 
5 representation of x including a binary exponent e and the binary mantissa m, wherein a 
first, most significant part of the partition corresponds to a region i and a second, less 
significant part of the partition corresponds to a region Ax, where Ax is a distance 
, i + 0.5 

from mantissa m to reference point a, = H — — ; and 

N 

compute an approximation to log(x), using a polynomial of first degree 
10 in m and a precomputed value of log(«j). 

25. A computing device in accordance with Claim 24 wherein said 
device being configured to compute the approximation to log(x) comprises said device 
being configured to compute an approximation written as: 

log(m) « log(a ; ) + — — — ; 

a i 

1 5 where a t is a closest reference point to the mantissa m; and 

l<a ; <2. 

26. A computing device in accordance with Claim 24 wherein said 
device being configured to compute an approximation to log(x) comprises said device 
being configured to compute an approximation written as: 

20 y = - log(x) » b l + c, Ax + e x log(2) 

for / = 0,...,iV-l 

where: 



-15- 



15-CT-5271 



b. =-iog(a,) + (— — I -[ i + -tt|— ; and 

&K ,J {4 ai N J { 2NJa i 
Ci =-\la ; . 

27. A computing device in accordance with Claim 26 further 
configured to precompute a value for log(2), and, for each i, to precompute each value 
of bi and c ; . 

28. A computing device in accordance with Claim 27 further 
configured to store the precomputed values of b { and c t in a look-up table. 
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METHOD AND APPARATUS FOR FAST 
NATURAL LOG(X) CALCULATION 

ABSTRACT OF THE DISCLOSURE 

The present invention is, in one embodiment, a method for computing 
a natural logarithm function that includes steps of: partitioning a mantissa region 
between 1 and 2 into N equally spaced sub-regions; precomputing centerpoints a t , of 
each of the N equally spaced sub-regions, where i=0,...,N-l; selecting N sufficiently 
5 large so that, within each sub-region, a first degree polynomial in m computes log(m) 
to within a preselected degree of accuracy for any m within the sub-region, where m is 
a mantissa of a binary floating point representation of a number; and computing a 
value of log(x) for a binary floating point representation of a particular number x 
stored in a memory of a computing device utilizing the first degree polynomial in m. 

10 This embodiment of the present invention and others described herein 

reduce the complexity of approximations used to calculate natural logarithms while 
achieving numerical accuracy consistent with IEEE floating point precision. 
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as stated below next to my n 



As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship a 

I believe I am the original, first and sole inventor (if only one name is listed below) or an original, first and joint inventor 
(if plural names are listed below) of the subject matter which is claimed and for which a patent is sought on the 
invention entitled: METHOD AND APPARATUS FOR FAST NATURAL 10G(X) CALCULATION, the specification of which: 



s attached hereto 
was filed on 



n Serial No. 



I hereby state that I have reviewed and understand thf 
claims, as amended by any amendment referred to above. 



I acknowledge the duty to disclose information which i 
with Title 37, Code of Federal Regulations §1 .56(a). 



of the above identified specification, including the 
material to the examination of this application in accordance 



1 hereby claim priority benefits under Title 35, United States Code, §120 of any United States application(s) listed below 
and insofar as the subject matter of each of the claims of this application is not disclosed in the prior United States 
application in the manner provided by the first paragraph of Title 35, United States Code, §112. I acknowledge the duty 
to disclose material information as defined in Title 37, Code of Federal Regulations, § 1.56(a) which occurred between 
the filing date of the prior application and the national or PCT international filing date of this application: 



Application Serial No. 



Filing Date 



Status (patented , pending, abandoned) 



I hereby claim the benefit under Title 35, United States Code §1 19(e) of any United States provisional application(s) listed 
below: 

Application Serial "No. Filing Date Additional provisional application 

numbers are listed on a 
" — supplemental 

Z^H^^^ZZZ^ZZZ^ priority sheet attached hereto. 



POWER OF ATTORNEY: As a named inventor, 1 hereby appoint the following attorney(s) and/or agent(s) to prosecute 
this application and transact all business in the Patent and Trademark Office connected therewith, (.list name and 
registration number) 

John S. Beulick (Reg. No. 33,338), and Alan L. Cassel (Reg. No. 35,842), of Armstrong 
Teasdale LLP, One Metropolitan Square, Suite 2600, St. Louis, MO 63102-2740; Phyllis Y. 
Price (Reg. No. 34,234) and Christian G. Cabou (Reg. No. 35,467), of GE Medical Systems, 
3000 North Grandview Blvd., W-710, Waukesha, Wl 53188; Jay Chaskin (Reg. No. 24,030); 
Ronald E. Myrick (Reg. No. 26,315) and Henry J. Policinski (Reg. No. 26,621), of General 
Electric Company (W3C), 3135 Easton Turnpike, Fairfield, CT 06431-0001 


Send Correspondence to: 


Direct Telephone Calls To: 


John S. Beulick 


John S. Beulick 


Armstrong Teasdale LLP 


314/621-5070 


One Metropolitan Square, Suite 2600 




St. Louis, MO 63102-2740 
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I hereby declare that all statements made herein of my own knowledge are true and that all statements made on 
information and belief are believed to be true; and further that these statements were made with the knowledge that 
willful false statements and the like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful false statements may jeopardize the validity of the application 
and any patent issued thereon. 



SOLE OR FIRST INVENTOR: 



Signature: Date: 

Residence: Waukesha, Wisconsin 53189 

Citizenship: _ P.R. China 

Post Office Address: _ 2213 Stony Ridge Drive, Waukesha, Wisconsin 53189 



SECOND JOINT INVENTOR, IF ANY: 

Full Name: Sharon X. Wang 

Signature: Date: 

Residence: Brookfield, Wisconsin 53045 

Citizenship: P.R. China 



Post Office Address: 3325 Mountain Drive, Brookfield, Wisconsin 53045 



